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Abstract 

Assuming positive observation of neutrinoless double bata decay together 
with the CHOOZ reactor bound, we derive constraints imposed on neutrino 
mixing parameters, the solar mixing angle 6\2 and the observable mass param- 
eter (m)/3 in single beta decay experiments. We show that 0.05 eV < {m)p < 
2 eV at the best fit parameters of the LMA MSW solar neutrino solution by 
requiring the range of the parameter (m)^ deduced from recently announced 
double beta decay events at 95 % CL with ±50 % uncertainty of nuclear 
matrix elements. 
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I. INTRODUCTION 



While the good amount of evidences for the neutrino mass and the lepton flavor mixing 
have been accumulated |l|-f|, we still lack observational indications of how large is the 
absolute mass of the neutrinos. To our understanding to date it may show up in only a few 
places, the single j| or the double beta decay || experiments as well as future cosmological 
observations ||. Other potential possibilities for hints of absolute mass of neutrinos include 
Z-burst interpretation of highest energy cosmic rays . 

Among these various experimental possibilities the neutrinoless double beta decay exper- 
iments seems to have relatively higher sensitivities. The most stringent bound on effective 
mass parameter {m)pp (see eq. (0) for definition) is now {m)pp < 0.35 eV, which comes from 
Heidelberg- Moscow group [EJ. Furthermore, a wide variety of proposals for future facilities 
as well as ongoing attempt with greater sensitivies are actively discussed. They include 
NEMO 0, GENIUS |T§, CUORE |TT|, MOON 0, XMASS |T3], and EXO projects. 
These high-sensitivety experiments open the enlighting possibility of discovering neutrinoless 
double beta decay events, not just placing an upper bound on (m)pp by its nonobservation. 
Therefore, it is of great importance to completely understand what kind of informations can 
be extracted if such discovery is made. 

We discuss in this paper in a generic three flavor mixing framework the constraints on 
neutrino masses and mixing by positive observation (as well as nonobservation) of neu- 
trinoless double beta decay. The constraints imposed on neutrino mixing parameters by 
neutrinoless double beta decay have been discussed by many authors. They include the 
ones in early epoch jnj, those in "modern era" in which real constraints on solar mixing 
parameters are started to be extracted |T^JT7H , and the ones in "post-modern era" where the 
analyses are performed in a complehensive manner in the framework of generic three flavor 
neutrino mixing |[18|| . 

In a previous paper, we have made a final step in the series of analyses by proposing a 
way of expressing the constraints solely in terms of observables in single and double beta 
decay [|H|. By using the framework, we discussed the possibility of placing lower bound on 
\U e ^\ 2 assuming positive observation in direct mass measurement in single beta decay and 
an upper limit on (m)pp in double beta decay experiments. It is a natural and logical step 
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for us to examine next the alternative case of positive observation of neutrinoless double 
beta decay events. 

Timely enough, an evidence of the neutrinoless double beta decay has just been reported 
by Klapdor-Kleingrothaus and collaborators [ 20| . Since the confidence levels of the claimed 
evidence are about 2 and 3 a in Bayesian and Partcle Data Group methods, respectively, we 
must wait for confirmation by further data taking, or by other groups to conclude that neu- 
trinos are Majorana particles. Nevertheless, we feel that the peak in the relevant kinematic 
region in their experiments is too prominent to be simply ignored. 

As will become clear as we proceed it is essential to combine the constraint on |?7 e 3| 2 = s^ 3 
imposed by the reactor experiments |^TJ . One of the key points in our subsequent discussion 
is that the double beta and the reactor bounds cooperate to produce a stringent constraint 
on absolute mass scale of neutrinos and the mixing angle 9 12 which is responsible for the 
solar neutrino problem. 



II. CONSTRAINTS FROM NEUTRINOLESS DOUBLE BETA DECAY 



Let us start by defining our notations. We use throughout this paper the standard 



notation of the MNS matrix [22 



U 



C12C13 



S12C13 



si 3 e 



-S12C23 - Ci2S 2 3Sl3e 4(5 C12C23 - Si2S23Si 3 e i<5 S 23Cl3 



S12S23 - Ci 2 C23Si 3 e 



-C12S23 - si 2 C2 3 si 3 e C23C13 



(1) 



Using the notation, the observable in neutrinoless double beta decay experiments can be 
expressed as 



(m) 



i=l 



(2) 



where m; (i=l, 2, 3) denote neutrino mass eigenvalues, U e i are the elements in the first 
low of the MNS matrix, and /3 and 7 are the extra CP-violating phases characteristic to 
Majorana neutrinos ||23|| , for which we use the convention of Ref. (JTB] . 
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We define the neutrino mass-squared difference as Am?- = m 2 — mf in this paper. In 
the following analysis, we must distinguish the two different neutrino mass patterns, the 
normal (Am| 3 > 0) vs. inverted (Am| 3 < 0) mass hierarchies. We use the convention that 
W3 is the largest (smallest) mass in the normal (inverted) mass hierarchy so that the angles 
9 12 and # 23 are always responsible for the solar and the atmospheric neutrino oscillations, 
respectively. We therefore often use the notations |Am| 3 | = Am 2 trn and Am^ = Am|, to 
emphasize that they are experimentally measurable quantities. Because of the hierarchy of 
mass scales, Amg/Am^ m <C 1, Am^ can be made always positive as far as 9 12 is taken in 
its full range [0, tt/2] f2|. 

In order to derive constraint on mixing parameters we need the classification. 

> m 3 s\ 3 (3) 

< ™3^ 3 (4) 

However, examination of the case II reveals that it does not lead to useful bounds. Therefore, 
we only discuss the case I in the rest of this paper. 



Case I: 



m i c i2 c i3 e 1(3 + m 2 s? 2 c? 3 e' 



-i/3 



Case II: 



m l C 12 C 13 e iP + m 2S 2 i2C 2 n e +iP 



A. Joint constraint by upper bounds on {m)pp and reactor experiments 

Since we try to utilize the experimental upper bound on (rn)pp, (m)pp < ( m ) r pp x \ we 
derive the lower bound on (m)pp. It can be obtained in the following way; 



(m)pp > c 2 13 



13 



'm\(? X2 + m 2 s i2) cos P ~ i(fnicl 2 ~ m 2sl 2 ) sm P 
m i c i2 + m 2 s i2 + 2mim 2 cf 2 Si2 cos 2(3 — m 3 sf 3 . 



m 3 s\ 3 



(5) 



Noticing that the right-hand-side (RHS) of (|5|) has a minimum at cos 2/5 = —1, we obtain 
the inequality 



(m)pp > c 



13 



m\(? 12 — rri2s\ 2 



m 3 s 



13- 



(6) 



We note that the RHS of (||) is a decreasing function of s 2 3 , and hence takes a minimum 
value for the maximum value of s^ 3 which is allowed by the limit placed by the reactor exper- 



iments [ 21 1 . We denote the maximum value as s^ 3 (CH) throughout this paper. Numerically, 



sf 3 (CH) ~ 0.03. (While the precise value of the CHOOZ constraint actually depends upon 
the value of Am^ m |2T] , we do not elaborate this point in this paper.) Using the constraint 
we obtain 



(m) 



m\(? 12 — m 2 s\ 2 



- (rn 3 + 



m\(? 12 — m 2 s\ 2 



4(CH). 



(7) 



It can be rewritten as the bound on cos 2#i 2 = cos 29 Q as 

m 2 - mi (m)^ + m 3g f 3 (CH) ^ m 2 - m, jmgHn^Cg 

T7 \ 9 /riTTN — COS ^12 S 1 T7 \ 9 /,~itt\ 

m 2 + m 1 i(m 2 + mi)cf 3 (CH) m 2 + m 1 ±{m 2 + mi)cf 3 (CH) 

where c? 3 (CH) = 1 -s? 3 (CH). 



(8) 



B. Joint constraint by lower bounds on (m)^ and reactor experiments 

A positive observation of neutrinoless double beta decay will lead to the experimental 
lower bound on (m)pp, > ( m )/?/3™> which we use to place new bound on neutrino 

mixing parameters. Toward the goal we note, similarly as (|Sp, that 

{m)pi3 < c i3\/ m i c i2 + m hi2 + 2mim 2 Ci 2 sj 2 cos 2(3 + m 3 sl 3 , (9) 

whose RHS is maximized by taking cos 2(3 = +1 and sf 3 = s^ 3 (CH) in the last term and 
cf 3 = 1 in front of the square root. (A more refined treatment entails the same excluded 
region.) One can then derive an inequality similar to (0); 

(m)™" < (m)p < (m x c? 2 + m 2 s\ 2 ) + m 3 ^ 3 (CH). (10) 

By rewriting (pT0|) we obtain the other upper bound on 00826*12; 

cos26> 12 < ^7 r . (11) 

m 2 -mi ±(m 2 -m 1 ) 

To summarize, we have derived in this section the two kinds of upper bound on cos 2# 12 
(lower bound for cos2# 12 < 0) by using the assumed experimental constraint (m)™™ < 
( m )i3i3 < (m)™p x imposed by neutrinoless double beta decay experiments. 
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III. CONSTRAINTS EXPRESSED BY EXPERIMENTAL OBSERVABLES 



We rewrite the bounds on solar mixing angle in terms of measurable quantities. Toward 
the goal we note that three neutrino masses rrii (i=l,2,3) can be expressed by the two 
Am 2 and a remaining over-all scale m#. We assign m# to the mass of the highest-mass 
state, 777,3 in the normal mass hierarchy (Am^ > 0), and m<i in the inverted mass hierarchy 
(Am 2 z < 0), respectively. 



We have argued in our previous paper JT9| that in a reasonable approximation one can 
regard m# as the observable (m) p in direct mass measurements in single beta decay experi- 
ments!] We have noticed that the identification is exact in two extreme cases of degenerate 
and hierarchical mass spectra. Then, the three mass eigenvalues of neutrinos can be repre- 
sented solely by observables; Am^ im , Arrig, and (m)p in a good approximation. 

In each neutrino mass pattern, we have the expressions of three mass eigenvalues: 
Normal mass hierarchy (Am| 3 > 0); 



m i = y m H ~ Am 2 atm - Ami, m 2 = V m 2 H - Am 2 atm , m 3 = m H ~ (m) p. (13) 
Inverted mass hierarchy (Am^ < 0); 



m i = V m H ~ Am !> m 2 = m H ~(m) /3 , m 3 = \Jm 2 H - Am 2 atm . (14) 

It is instructive to work out the form of constraint in the degenerate mass approximation, 
m 2 ~ m 2 ^> Am 2 tm , Am 2 Q . It is easy to show that in the degenerate mass limit the bound 
(|) becomes 



cos2fl 12 | < sec 2 fl 13 (CH) 



+ 4(CH) 



(15) 



( m )l3 

On the other hand, the bound @ gives the inequality (m)p > (m)™™ in the degenerate 
mass limit. (To show this one may go back to (|^), rather than using (|TT|).) 



*While we used the linear formula derived by Farzan, Peres and Smirnov |2 



<-j i 



V n \TT -I 2 
2^j=l \ u ej\ 



(12) 



with n being the dimension of the subspace of (approximately) degenerate mass neutrinos, this 



point remains valid even if we use an alternative quadratic expression 26]. 



IV. ANALYSIS OF THE DOUBLE BETA-REACTOR JOINT CONSTRAINTS 



We analyze in this section the joint constraints derived in the foregoing sections and 
try to extract the implications. Let us start by examining the case of recent obsevation 



announced in [2(J which gives rise to 0.11 eV < (m)pp < 0.56 eV and 0.05 eV < (rn)pp < 
0.84 eV if ±50 % uncertainty of the nuclear matrix elements are considered, each at 95 % 
CL. In Fig. 1 we present on (m)p - cos2# 12 plane the constraint (M) by the thick solid 



lines (solid line) and (|TT|) by the thick dashed line (dashed line) for cases with (without) 
uncertainty of the nuclear matrix elements, respectively. The regions surrounded by these 
lines are allowed. The slope of (m) i g-dependence of ( p|) is so large that the dashed line 
looks like a vertical line, which implies the inequality (m)p > (m)p n p n . We have derived it 
earlier in the degenerate mass limit, but it is generically true if AttIq is smaller than other 
relevant mass squared scales. Only the case of normal mass hierarchy (Am^ > 0) is shown 
in Fig. 1; the case of inverted hierarchy (Am| 3 < 0) gives an almost identical result except 
for a slight shift of the dashed line toward smaller (m)p by ~ 10 %. 

Superimposed in Fig. 1 are the 95 % CL allowed regions of cos 2# 12 for the large mixing 
angle (LMA) MSW solution (indicated by the shaded region between thin solid lines) and 
the low (LOW) MSW solution (indicated by the shaded region between thin dashed lines) 
of the solar neutrino problem |?7]] . There are several up to date global analyses of the solar 
neutrino data with similar results of allowed region of mixing parameters |p8|j . Therefore, 



we just quote the result obtaind by Krastev and Smirnov in the last reference in |28[ . 

Figure 1 illustrates that for a given value of cos 2# 12 the single beta decay observable (m) p 
has to fall into a region bounded by (m)™ m ~ (m)™™ and (m)g ia:E , which are dictated by ( |TTD 
and (|D, respectively. Thus, we have a clear prediction for direct mass measurements using 
a single beta decay with observation of double beta decay events. With use of the numbers 
given in ||20|| , for example, the observable {m)p must fall into the region 0.05 eV < (m)p < 
2 eV (0.11 eV < (m)p < 1.3 eV) with (without) uncertainty of nuclear matrix elements at 
the best fit parameters of the LMA MSW solution. (The best fit value is tan 2 (9 12 = 0.35, 
or cos2# 12 = 0.48 in the last reference in ||28|| .) Within the allowed region the cancellation 
between three mass eigenstates can take place for appropriate values of Majorana phases 
that allow (typically) a factor of 2-3 larger values of (m)p compared with the measured value 
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of {m)pp. At around maximal mixing (cos2#i2 — 0), which is allowed by 95 % CL in the 
LOW solution, the cancellation is so efficient that much larger values of (m)p is allowed. 
Therefore, there are still ample room for hot dark matter mass neutrinos both in the LMA 
and the LOW solutions. 

It should be emphasized that a finite value of {w)pp does imply a lower bound on (m)p, 
as indicated in Fig. 1; a vanishingly small (m)p cannot be consistent with finite (m)pp in 
double beta decay experiments. The sensitivity of the proposed KATRIN experiment is 
expected to extend to (m)p < 0.3 eV P9| . On the other hand, the present 68 % CL limit 
quoted in |20| without nuclear element uncertainty is 0.28 eV < {m)pp < 0.49 eV. Therefore, 
if the limit is further tightened by additional data taking in the future, both experiments 
can become inconsistent, giving an another opportunity of cross checking. 

In Fig. 2, we demonstrate the approximate scaling relation obeyed by the constraint @ 
by taking (m)p/(m)^ x as the abscissa in a wide range of the (m)p n p x in degerarate mass 
region, 0.1 eV < (m)^ x < 1 eV. The scaling is exact in the degenerate mass limit as shown 
in (IT5J). The relation is useful to estimate the allowed region of (m)p for a given value of 
(m)™p x which is not explicitly examined in this paper. 

The most stringent bound to date on {m)p is from the Mainz collaboration |p(J|| , ( m )i3 < 
2.2 eV (95 % CL). (A similar bound (m)p < 2.5 eV (95 % CL) is derived by the Troitsk group 
f3~T| .) As we can see in Fig. 2 that the double beta bound with the CUORE sensitivity region 
(m)pp <0.3 eV Jill becomes stronger than the Mainz bound for the LMA MSW solution 
but not for the LOW MSW solution in their 95 % CL regions. 

In Fig. 3, we present the similar allowed regions for hypothetical discovery of neutrinoless 
double beta decay events which would produce the experimental bounds 0.01 eV < (m)pp < 
0.03 eV. It is to examine how the constraint changes in some other situation of discovery 
with different mass parameter ranges. We note that even such deep region of sensitivity will 
be explored by several experiments PH|JT^ -P . 



For this case, the bounds for the normal and the inverted mass hierarchies start to split as 
shown in Fig. 3. In the case of inverted mass hierarchy the lower bound on {m)p is replaced 



by the trivial bound (m)p > J /\m\ tm which is more restrictive. The latter is indicated by 
the dash-dotted line in Fig. 3b. It is also evident that the constraint from double beta decay 



S 



is so stringent that the limit on (m)^ is tightened to be (rn)p< 0.2 eV for the LMA MSW 
solution. 

In conclusion, we have demonstrated in this and the previous papers the mutual intimate 
relationship between observation and/or nonobservation in single beta decay and neutrino- 
less double beta decay experiments. We hope that it stimulates even richer future prospects 
not only in double beta decay experiments but also in direct mass measurements using single 
beta decay. 

Finally, some remarks are in order: 



(1) If the LMA MSW solution is the case and if the KamLAND experiment [32) that just 
started data taking can measure cos26 , 12 within 10 % level accuracy, the upper limit of (m)^ 
can be accurately determined with ~ 20 % accuracy. 

(2) In this paper we have derived constraints imposed on neutrino mixing parameters by 
observation of neutrinoless double beta decay events and the CHOOZ reactor bound on 
| Z7e3 1 2 in the generic three flavor mixing framework of neutrinos. Suppose that neutrinoless 
double beta decay events are confirmed to exist and the single beta decay experiments detect 
neutrino mass outside the region of the bound derived in this paper. What does it mean? 
It means either that double beta decay would be mediated by some mechanisms different 
from the usual one with Majorana neutrinos, or the three flavor mixing framework used in 
this paper is too tight. 

Nore added: 

After submitting the first version of our paper to the electronic archive, we became aware 
of the works which address relatively model-independent implication of the results reported 
in ||20|| , or critically comment on the interpretation of the events. References [33j and [34 
are the incomplete lists of them. 
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FIG. 1. The constraints imposed on neutrino mixing parameters 612 and the observable mass 
parameter {m)p in single beta decay experiments by recent observation of neutrinoless double bata 
decay events. The solid and the dashed lines represent the bounds @ and (|l~l|), respectively; the 
allowed region is inside these three lines. The bold and the normal lines are for the ranges of mass 
parameter 0.05 eV < {m)pp < 0.84 eV and 0.11 eV < {m)pp < 0.56 eV corresponding, respectively, 
with and without ±50 % uncertainty of nuclear matrix elements. The mixing parameters are fixed 
as Ara^ (m = 3 x 10 -3 eV 2 and Attiq = 4.8 x 10~ 5 eV 2 . Also shown as shaded region are the 
allowed regions of cos 2#i2 at 95 % CL for the LMA (the region between thin solid lines) and LOW 
(the region between thin dashed lines) MSW solutions. 
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FIG. 2. The approximate scaling relation obeyed by the bound (|8|) for a wide range of {m)^ x , 
0.1 eV (dashed lines) < (m)™p x < 1 eV (solid lines). The bound (|ll|) is schematically drawn by 
vertical thin dashed line. The shaded regions are the same as in Fig. 1. 
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FIG. 3. The same as in Fig. 2 but with assumed observed mass parameter (m)pp in the range 
0.01 eV < (m)j3/3 < 0.03 eV. Fig. 3a and 3b for the normal and the inverted mass hierarchies, 
respectively. The dash-dotted line in Fig. 3b denotes the trivial bound (m)^ > J Am„ (m . 
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